Spectral functions in itinerant electron systems with geometrical frustration 
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The Hubbard model with geometrical frustration is investigated in a metallic phase close to half- 
filling. We calculate the single particle spectral function for the triangular lattice within dynamical 
cluster approximation, which is further combined with non-crossing approximation and fluctuation 
exchange approximation to treat the resulting cluster Anderson model. It is shown that frustra- 
tion due to non-local correlations suppresses short-range antiferromagnetic fluctuations and thereby 
assists the formation of heavy quasi-particles near half-filling. 



Recently, geometrically frustrated metallic systems 
have attracted much attention, for which a wide variety 
of interesting phenomena have been discovered. For in- 
stance, the compound L1V2O4 with pyrochlore structure, 
which is given by a corner-sharing network of tetrhedra 
of V ions, exhibits a heavy fermion behaviora. Also, an- 
other pyrochlore compound Y(Sc)Mn20 shows a quantum 
spin-liquid behavior down to low temperatures. These in- 
teresting phenomena are cnHgideced to be closely related 
to geometrical frustratiorJj'QtfHiJ'El induced by the strong 
Coulomb interaction with specific geometry of the lattice. 

Among theoretical approaches to strongly correlated 
electron systems, dynamical mean field theory (DMFT)_- 
is one of the most successful methods to describe various 
physical properties. In DMFT the self energy is given as 
a local quantity, which is justified in the limit of large 
dimensionsul__. This method has been known as a good 
approximation even in three dimensions. However, since 
the self energy becomes local in DMFT, non-local charge 
and spin correlations, which are essential for treating geo- 
metrical frustration, cannot be properly described in this 
framework. In order to deal with such non-local cor- 
relationS|-4jafiipu8-|a«proaches beyond DMFT have been 
proposecOll3'til'E3ll3. _A.mpn.ffi others, dynamical cluster 
approximation (DCA)0_3i__'__O may provide a system- 
atic way to incorporate non-local correlations. In DCA, 
the lattice problem is replaced by the corresponding clus- 
ter one embedded in an effective medium determined self- 
consistently. This method has the following nice features: 
the algorithm is fully causal and the approximation can 
be improved systematically if the cluster size is increased. 

In this paper, we investigate the effects of geometri- 
cal frustration on the Hubbard model. In particular, we 
focus on the triangular lattice, which is known as a pro- 
totypical system having strong geometrical frustration. 
Since it is important to incorporate short-range correla- 
tions systematically, we employ DCA, which is fur-thpr. 
combined with non-crossing approximation (NCA)I__'E_I 
to solve the local problem. We also use a weak-coupling 
approach byjaeans of fluctuation exchange approxima- 
tion (FLEX)I__ for the local problem, which should be 
complementary to the treatment of NCA. By calculating 
the one-particle spectral function as well as the total den- 
sity of states in these approximations, we discuss how ge- 
ometrical frustration affects the formation of heavy quasi- 
particles in a metallic phase close to half-filling. 



We start with the single-band Hubbard model, 

H = -tJ2 chcjv-t' 4a c f'r + U^2n n n ii , (1) 
{i,i)<* {*•?)<' * 

where Cj-(cj ) is the annihilation (creation) operator of 
an electron with spin a at the i-th site, and U repre- 
sents the Coulomb repulsion. We introduce two kinds of 
hopping parameters t and t'(< t) to study the effects of 
geometrical frustration svstematicallvEl (FIG. |l| (a)). 




FIG. 1: (a) Schematic representation of triangular lattice with 
electron hoppings t and t' . (b) Example of the coarse-graining 
cells in the BZ for the triangular lattice, where the cluster 
size N c —4 and the dashed line denotes the first BZ. The dots 
represent the cluster momenta K. 

In order to apply DCA to the frustrated lattice sys- 
te m ,) mc ijbg^introduce the N c discrete cluster-momenta 
which are defined as, 



K = rtai + m&2 
ai = 



I, 83= 0, 
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(2) 
(3) 



where n and m denote the integral number and K should 
be within the first Brillouin zone (BZ). Next, we di- 
vide the original BZ into the subregions specified by 
the cluster momenta (coarse-graining cells); an exam- 
ple of N c = 4 is shown schematically in Fig. |l| (b). 
Therefore., in pu| treatment the coarse-grained Green 
function! 17 ''' 
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is specified by the cluster momentum K. Here, N' — 
N/N c where N is the number of total lattice sites. The 
summation over k is taken within the coarse-graining cell. 

In order to obtain the coarse-grained Green function, 
we now map th@-,Hubbard lattice model to the cluster 
Anderson modefij. By solving this effective cluster prob- 
lem, we can obtain the cluster self energy E(K, z), so 
that the coarse-grained Green function is determined self- 
consistently. Then the one-particle spectral function and 
the total DOS are given by the standard formula, 

A a (K,uj) = -ilmG CT (K» (5) 

c K 

We numerically iterate the above procedure until the 
calculated quantities converge within desired accuracy. 
In the following discussions, we shall deal with a param- 
agnetic metallic phase. We set t = 1 as the unit of the 
energy for simplicity. 

In order to solve th^cluster problem mentioned above, 
we make use of NCAEj, which is expected to provide re- 
liable results in the temperature range we are now inter- 
ested in. Since the dimension of the cluster Hamiltonian 
is 4 Na , so that the numerical calculation becomes much 
more difficult with the increase of the number of cluster 
momenta. We practically take the cluster size, N c = 4 
within NCA in this paper. Note that we have neglected 
off-diagonal terms in the cluster resolvents, since those 
are known to b£_less important for discussing the one- 
particle spectraO. We show our results calculated for 
T = 0.6 below, since this temperature is reasonably low 
for our system to exhibit essential properties due to the 
electron correlations. 

Let us first discuss the total DOS in the system with 
strong frustration. In FIG. @, we show the DOS for the 
triangular lattice at half-filling with isotropic hopping, 
t' = t. For the triangular lattice, DOS calculated by 
DMFT has a many-body peak around the Fermi level 
(u> ~ 0), which implies the formation of heavy quasi- 
particles due to the Hubbard interaction. We note here 
that the DC A recovers the DMFT for N c = 1. Even if 
short-range non-local correlations are taken into account 
by DCA, the many-body peak still persists although it 
becomes slightly smaller than that of DMFT. 

This Fermi-liquid like behavior is contrasted to the 
case of the square lattice with nearest neighbor hopping. 
In the latter case, it is known that a pseudo-gap struc- 
ture in DOS, which is developed by the Hubbard- inter- 
action, prevents the formation of quasi-particlest^l when 
the system is at half-filling where short-range antiferro- 
magnetic fluctuations are most relevant (see also FIG. 
^). Our results in the triangular lattice shows that anti- 
ferromagnetic correlations are strongly suppressed in the 
frustrated lattices, resulting in the formation of heavy 
quasi-particles even at half filling. Note that in both of 
these two lattices, the Mott insulator is stabilized as U 
becomes large. 




FIG. 2: Total DOS for the triangular lattice with DCA (solid 
line) and DMFT (dashed line) at half-filling, respectively. 
Here, the cluster size N c = 4. The energy is measured from 
the Fermi level (u = 0). The parameters are chosen as t' = t, 
(7=6.0 and the temperature T=0.6 (t is taken as the unit of 
the energy). The inset shows DOS for the non-interacting 
tight-binding model on the same lattice. 



To investigate how antiferromagnetic fluctuations are 
suppressed by geometrical frustration, we change t' con- 
tinuously, and observe what happens for DOS and the 
one-particle spectral functions on the triangular lattice. 
The results are shown in FIG. Note that the effect of 
frustration becomes less important with the decrease of 
t' , As t' decreases from t, the heavy quasi-particle band 
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FIG. 3: (a) Total DOS p(ui) and (b)-(d) one-particle spectral 
functions A (K, w) corresponding to the K = (0,0), (ir, 7t/v3) 
(or (tt, -tt/V3)) and (0, 2tt/v / 3) for various choices of the hop- 
ping amplitude: t' — t (solid line), t' — 0.5t (dashed) and 
t' — 0.0 (dotted). The other parameters are as in FIG. H. 
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around the Fermi level is obscured, and then the pseudo- 
gap is developed after frustration is suppressed. From 
the one-particle spectral function for each momentum, 
shown in FIG. || (b)-(d), we can figure out which con- 
tribution is most relevant to the quasi-particle formation 
(or suppression) under strong frustration. Since the state 
with K = (0, 0) is away from the Fermi level, the overall 
structure is not changed by t' . However the spectrum for 
K = (tt, 7t/v3), which is energetically degenerate with 
that for K = (tt, —k/VS), crosses the Fermi surface, be- 
ing easily affected by frustration. When t' is equal to 
t, where frustration is very strong, heavy quasi-particles 
are formed around the Fermi level. With decreasing frus- 
tration, the heavy quasi-particle band splits, and the 
pseudo-gap begins to develop explicitly around the Fermi 
level. Therefore, the low energy physics in this system is 
mainly controlled by the state with K = (tt, n/^/S) (or 
K=(7r,-7r/V3)). " 

In order to further confirm the above statements, we 
have numerically estimated the renormalization factor, 
which is defined as, 



Z K = | 1 - 



95RS(K,w) 



du> 



(7) 




FIG. 4: Renormalization factor for the triangular lattice, (a) 
[/-dependence (t! — t) and (b) t'-dependence (U = 6.0). Here, 
K = (tt,tt /\f3). The other parameters are as in FIG. [| 



anomalous behavior in Zk is closely related to the for- 
mation of a dip-structure in the spectrum, which reflects 
the enhancement of antiferromagnetic fluctuations. The 
negative values of Zk should be considered as an artifact 
of fitting the numerical data with the formula (Q) even 
for the case with a dip structure. 

So far, we have treated the effective cluster model by 
NCA, which may be an efficient approximation in the 
intermediate coupling regime. Here, we take a weak- 
coupling approach using FLEX to solve the effective clus- 
ter problem. Since FLEX is a perturbative method with 
respect to the Coulomb interaction U, it may give re- 
sults complimentary to those of NCA. We here confirm 
the conclusion of the NCA approach by investigating ef- 
fects of frustration on the triangular lattice in Fig.[j](a). 
Shown in Fig.|| are the total DOS, the one-particle spec- 
tral function and the corresponding self-energy calculated 
for the triangular lattice at half-filling. 
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In FIG. |] we show the results obtained for the mo- 
mentum K = (tt, tt/^/S) close to the Fermi level. It 
is seen from FIG. [|(a) that the renormalization factor 
for the triangular lattice is reduced with the increase of 
U for the fully frustrated case (f — t). As mentioned 
above, this indicates that the effective mass near the 
Fermi level is enhanced by U, forming well-defined heavy 
quasi-particles. When t' /t is decreased from unity with U 
being fixed, the effect of frustration becomes less impor- 
tant, so that the system has a tendency to be an insulator 
with strong antiferromagnetic fluctuations. In this case, 
we encounter an anomalous behavior in the renormal- 
ization factor calculated numerically for the triangular 
lattice, as seen from FIG.|j. Namely, the renormaliza- 
tion factor increases as t' decreases, and diverges around 
t'/t ~ 0.4, below which it takes negative values. This 



FIG. 5: (a) DOS p(ui), (b) one particle spectral function 
A(K.,uj) (K = (7r,7r/\/3)) and (c) imaginary part of the self- 
energy E(K,w) for U = 3.0, T = 0.6 at half-filling. The 
cluster size N c used for the calculation is 36. 

In the case of t' = 0, a pseudo-gap develops around the 
Fermi level in the total DOS as well as in the one-particle 
spectral function &$JU. increases, which is consistent with 
the QMC resultsElHEl and also with those obtained by 
NCA. It is seen that the pseudo-gap disappears with the 
increase of t' , similarly to the results of NCA. However, 
it is not clearly seen from this figure whether a heavy 
quasi-particle band is indeed formed, since the bump 
structure in the DOS above the Fermi level (t' = 0.5£ 
and l.Ot) is much affected by the van-Hove singularity. 
Nevertheless, we can see a tendency to the formation of 
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Fermi quasi-particles by observing the imaginary part of 
the self-energy in the large t' regime. (Fig.g(c)) Namely, 
when t' is absent, the imaginary part of the self-energy 
has a sharp peak structure around u> ~ 0, which is quite 
different from the Fermi-liquid behavior. However, with 
increasing t', its shape changes from the peak to a con- 
cave structure, which is similar to a conventional Fermi 
liquid-like behavior. We have also checked that qualita- 
tively analogous behaviors can be found in the system 
of different cluster size, N c — 16 and 64. Therefore, we 
confirm the conclusion of the NC A approach that frustra- 
tion assists the formation of Fermi quasi-particles by ef- 
fectively suppressing the antiferromagnetic fluctuations. 

We wish to mention here that our results are con- 
sistent with some experimental findings, e.g. the large 
specific-heat coefficient found in the pyrochlore com- 
pound, Y(Sc)Mn20, which is a typical example of geo- 
metrically frustrated metallic systems. 



In summary, we have investigated the Hubbard model 
with geometrical frustration by employing the triangu- 
lar lattice as a typical example. Applying DCA com- 
bined with NCA and FLEX to the Hubbard model in 
a paramagnetic metallic phase close to the Mott insu- 
lator, we have calculated the one-particle spectral func- 
tion. Based on the results obtained by both approaches, 
we have demonstrated how geometrical frustration sup- 
presses antiferromagnetic correlations, and then assists 
the formation of a heavy quasi-particle band near the 
Mott insulating phase. 

We acknowledge valuable discussions with A. Koga and 
A. Kawaguchi. The work is partly supported by a Grant- 
in- Aid from the Ministry of Education, Science, Sports, 
and Culture. Parts of the numerical computations were 
done by the supercomputer center at Institute for Solid 
State Physics, The University of Tokyo. 
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